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Free streaming in a mixture of collisionless non-relativistic dark matter (DM) particles is studied 
by solving the linearized Vlasov equation implementing methods from the theory of multicomponent 
plasmas. The mixture includes Fermionic, condensed and non-condensed Bosonic particles decou- 
pling in equilibrium while relativistic, heavy thermal relics that decoupled when non-relativistic 
(WIMPs), and sterile neutrinos that decouple out of equilibrium when they are relativistic. The dif- 
ferent components interact via the self-consistent gravitational potential that they source. The free- 
streaming length A/ 3 is obtained from the marginal zero of the gravitational polarization function, 
which separates short wavelength Landau-damped from long wavelength Jeans-unstable collective 

modes. At redshift z we find jrj^j = [t^t] 2 J2 a v °- 9d, a ("WkeV) 2 I a , where < v a < 1 are 

£SJ , the fractions of the respective DM components of mass m a that decouple when the effective number 

. of ultrarelativistic degrees of freedom is gd,a, and I a are dimensionless ratios of integrals of the dis- 

• tribution functions which only depend on the microphysics at decoupling and are obtained explicitly 

Qin all the cases considered. If sterile neutrinos produced either resonantly or non-resonantly that 
decouple near the QCD scale are the only DM component, we find A/ s (0) ~ 7kpc(keV/m) (non- 
resonant), A/ s (0) ~ 1.73 kpc (keV/m) (resonant). If WIMPs with m wimv > 100 GeV decoupling at 
t-H ' Td > lOMeV are present in the mixture with u wimp ^> 10~ 12 , A/ s (0) < 6.5 x 10 -3 pc is dominated 

by CDM. If a Bose Einstein condensate is a DM component its free streaming length is consistent 
' , ' with CDM because of the infrared enhancement of the distribution function. 

■ 
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Candidate dark matter (DM) particles are broadly characterized as cold, hot or warm depending on their veloc- 
, ity dispersions. The concordance ACDM standard cosmological model emerging from CMB, large scale structure 
J> ■ observations and simulations favors the hypothesis that DM is composed of primordial particles which are cold and 
collisionless [lj]. Although this model is very successful in describing the large scale distribution of galaxies, recent 
observations hint at possible discrepancies summarized as the "satellite" and "cuspy halo" problems. In the ACDM 
model the CDM power spectrum favors small scales, which become non-linear first and collapse in a hierarchical 
"bottom-up" manner and dense clumps survive the mergers in the form of "satellite" galaxies. Large-scale simula- 
tions within the ACDM paradigm lead to an overprediction of "satellite" galaxies |2j, which is almost an order of 
magnitude larger than the number of satellites that have been observed in Milky- Way sized galaxies 0, d, 0, H, @] • 

Furthermore, large-scale N-body simulations of CDM clustering predict a density profile monotonically increasing 
towards the center of the halos0, @, H, [§, U3|, with asymptotic behavior p(r) ~ r~ 7 where 1 < 7 < 1.5(1, H> G3 which 
describes accurately clusters of galaxies, but indicates a divergent cusp at the center of the halo. In contrast with 
these results, recent observations seem to indicate central cores in dwarf galaxiesfll], [IH, [TH, [3, [l5[, sparking the 
"cusps vs cores" controversy. 

Warm dark matter candidates (WDM) have been invoked as possible solutions to these potential discrepancies [l6l 
IT?}, these particles feature velocity dispersions intermediate between CDM and HDM and can relieve the satellite 
and cuspy halo problems. The clustering properties of collisionless DM in the linear regime depend on a fundamental 
length scale: the free-streaming length, that determines a cutoff in the power spectrum. Length scales larger than 
the free-streaming scale undergo gravitational instability, and shorter scales are damped. A simple estimate of the 
free-streaming length Xf s is obtained from the familiar Jeans' length by replacing the speed of sound by the velocity 
dispersion of the particle. An equivalent estimate is obtained by computing the distance that the particle travels 
within a dynamical (Hubble) time[l|[ A/ s ~ (V 2 )^ /H . 

However, a thorough assessment of DM particles and structure formation requires a more detailed and reliable 
determination of the free-streaming length. The necessity for this has been recently highlighted by the recent results 
of ref.[19j that suggest that the first stars form in filaments of the order of the free streaming scale. 
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Perturbations in a collisionless system of particles with gravitational interactions is fundamentally different from 
fluid perturbations in the presence of gravity. The (perfect) fluid equations correspond to the limit of vanishing mean 
free path. In a gravitating fluid pressure gradients tend to restore hydrostatic equilibrium with the speed of sound 
in the medium and short wavelength fluctuations are simple acoustic waves. For large wavelengths the propagation 
of pressure waves cannot halt gravitational collapse on a dynamical time scale. The dividing line is the Jeans length: 
perturbations with wavelengths shorter than this oscillate as sound waves, while perturbations with longer wavelength 
undergo gravitational collapse. 

In a gas of collisionless particles with gravitational interaction the situation is different since the mean free path is 
much larger than the size of the system (Hubble radius) and the fluid description is not valid. Instead the Boltzmann- 
Vlasov equation for the distribution function must be solved to extract the dynamics of perturbations 2(1 21 1. Just 



as in the case of plasma physics, the linearized Boltzmann-Vlasov equation describes collective excitations [22j. In the 
case of a collisionless gas with gravitational interactions these collective excitations describe particles free-streaming in 
and out of the gravitational potential wells of which they are the source. The damping of short wavelength collective 
excitations is akin to Landau damping in plasmas [22I] and is a result of the phenomenon of dephasing via phase 
mixing in which the particles are out of phase with the potential wells that they produce [23I]. This situation is similar 
to Landau damping in plasmas where dephasing between the charged particles and the self-consistent electric field 
that they produce lead to the collisionless damping of the collective modes [22|. For a thorough review of collective 
excitations and their Landau damping in gravitational systems see ref. [24| . 

Gilbert [251] studied the linearized Boltzmann-Vlasov equation in a matter dominated cosmology for non-relativistic 
particles described by an (unperturbed) Maxwell-Boltzmann distribution function. In this reference the linearized 
this equation was cast as a Volterra integral equation of the second kind which was solved numerically. The result of 
the integration reveals a limiting value of the wavevector below which perturbations are Landau damped and above 
which perturbations grow via gravitational instability. Eventually the redshift in the expanding cosmology makes 
wavevectors that are initially damped to enter the band of unstable modes and grow[25j. However, the dividing 
wavevector between damped and growing modes emerges at very early times during which the expansion can be 



neglected |25[. The results of the numerical study were consistent with replacing the speed of sound by the Maxwellian 
velocity dispersion in the Jeans length (up to a normalization factor of 0(1)). Gilbert's equations were used also 
by Bond and Szalay[26| in their pioneering study of collisionless damping of density fluctuations in an expanding 
cosmology. These authors focused primarily on massive neutrinos and solved numerically Gilbert's integral equation 
but approximated the Fermi-Dirac distribution function and provided a fitting function for the numerical results. 
Gilbert's equations were also solved numerically to study dissipationless clustering of neutrinos in ref. [27] but with 
a truncation of their Fermi-Dirac distribution function and an analytic fit to the numerical solution of the integral 
equation was provided. Bertschinger and Watts [2ll.[28| also studied numerically Gilbert's equation within the context 
of cosmological perturbations from cosmic strings and massive light neutrinos, and more recently similar integral 
equations were solved approximately for thermal neutrino relics in ref. [2^] . 

Most of the studies of the solutions of the linearized Boltzmann-Vlasov equation for collisionless particles addressed 
one single species of DM candidates 1 and generally in terms equilibrium distribution functions or approximations 
thereof, for example Maxwell-Boltzmann for relics that decouple non-relativistically or Fermi-Dirac (without chemical 
potential) (or truncations of this distribution) for neutrinos. 

However, it is likely that DM may be composed of several species, this possibility is suggested in ref. [13] and most 
extensions of the standard model generally allow several possible candidates, from massive weakly interacting particles 
(WIMPs) to "sterile" neutrinos ( SU{2) singlets) [H[ • 

Furthermore, several possible WDM candidates may decouple out of local thermodynamic equilibrium (LTE) with 
distribution functions which may be very different from the usual ones in LTE. This is the case for sterile neutrinos 
produced non-resonantly via the Dodelson-Widrow (DW) mechanism [32| or through a lepton-driven MSW (Mikheyev- 
Smirnov-Wolfenstein) resonance [33). 

A. Motivation and goals: In this article we study free streaming of decoupled collisionless non-relativistic (DM) 
candidates focusing on two aspects: 

• A mixture of (DM) particles including CDM and WDM candidates: typical studies of structure formation invoke 
either CDM or WDM, but it is likely that both candidates are present with different fractions v of the total 
DM component. In fact most particle physics extensions beyond the standard model have plenty of room for a 
variety of CDM, WDM or HDM candidates. Thus we allow all of these components, each one contributing an 



1 The exception being ref. |28H wherein massive neutrinos decoupled with a Fermi-Dirac distribution and a cosmic string source ex. <5 3 (r) 
were studied, neglecting any radiation component. 



3 



arbitrary fraction v to the total (DM) content of the Universe. Although the DM candidates are collisionless and 
do not interact directly with each other, they interact indirectly via the gravitational potential that they source. 
As a result the free-streaming length of the mixture is a non-trivial function of the individual free-streaming 
lengths. 

• Free streaming has mostly been studied in the above references with particles that decoupled either when 
ultrarelativistic (as is the case for neutrinos) or non-relativic as in the case of weakly interacting massive particles 
(WIMPs) but generally in local thermodynamic equilibrium (LTE), namely with Fermi-Dirac, Bose-Einstein or 
Maxwell-Boltzmann distributions respectively. We seek to obtain the corresponding free streaming lengths for 
particles that decoupled in or out of LTE with arbitrary isotropic distribution functions, without any truncation. 
This aspect is important for sterile neutrinos either produced non-resonantlv [32j or resonantlv[3l[ |33| because 
these particles decoupled while relativistic but out of LTE. Therefore, we consider the most general mixture of 
Fermionic and Bosonic thermal relics that decouple when relativistic, including the possibility of a Bose-Einstein 
condensate (BEC)[13|, heavy non- relativistic thermal relics, as the case of WIMPs, and sterile neutrinos that 
decoupled out of LTE when ultrarelativistic. 

In order to carry out this program analytically we first neglect the cosmological expansion and solve the linearized 
Boltzmann-Vlasov equation in the non-expanding case exactly by implementing methods from the theory of multi- 
component plasmas [22]. The neglect of the cosmological expansion is warranted by the detailed numerical study in 
refs. |25l. [26[ wherein it was found that the dividing wavevector between the Landau damped modes and the modes 
that grow under gravitational instability is insensitive to the expansion, just as in the case of the Jeans instability 
where the Jeans wavevector can be extracted in the non-expanding case include the redshift dependence of the density, 
speed of sound and scale factors a vosteriori [l8l] . The analytic exact solution for the free-streaming wave- vector today 
kf s (0) — 27r/A/ s (0) in terms of the full distribution functions without truncation, in or out of LTE is a main result 
of this program, one which yields a reliable determination of free-streaming lengths for mixtures of (DM) components 
that decoupled with arbitrary distribution functions. 

This program is carried out by implementing methods from the theory of multicomponent plasmas 22], in particular 
we obtain the "gravitational polarization" function]!!, [24| for a mixture of (DM) components akin to the dielectric 
response function of multicomponent plasmas [13] ■ The collective excitations are described by the zeroes of this function 
in the complex frequency plane and the free-streaming wave- vector kf s is identified as that wavevector that separates 
between the Landau-damped short wavelength modes and the gravitationally Jeans- unstable long-wavelength modes. 

Based on the Liouville evolution of the decoupled distribution functions and assuming that the expansion of the 
universe is slow enough so that it can be treated adiabatically, we provide a scaling argument that determines the 
following dependence of the free streaming length on the redshift, 

x fs (z) = x fs (o)VT+^. (l.i) 

B. Summary of Results: Our main result for the comoving free streaming length Xf s (z) at redshift z of mixed 
(DM) is 




where v a is the fraction of (DM) of each component with ^ Q v a = 1, a is the effective number of ultrarelativistic 
degrees of freedom at decoupling for each species (a) of mass m a , and the functions If,Ib are dimensionless ratios 
of integrals of the distribution functions of the decoupled particles which are determined by the microphysics at 
decoupling. Their explicit expressions in the cases considered are given in section (llllj) . The label F refer to all 
possible Fermions with chemical potential /i decoupled in LTE at a temperature Td while ultrarelativistic, and sterile 
neutrinos produced non-resonantly via the (DW) mechanism (32| for which the chemical potential vanishes and Lp[0] = 
21n(2)/3C(3) = 0.3844, and the label s refers solely to sterile neutrinos produced via a lepton-driven (MSW) resonance 
via the mechanism described in ref.[33]. The label B corresponds to condensed or non-condensed Bosons of mass 
m and chemical potential /i that decoupled at temperature Td while ultrarelativistic. The function Lb features an 
infrared divergence in the limits n/Td\ m/Td — » or fi = m for any value of the mass. This latter case corresponds 
to the case of a Bose-Einstein Condensate [34|. Thus Bosonic particles that decoupled while ultrarelativistic with 
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n/Td <C 1 or that formed a BEC lead to small free-streaming lengths and behave as CDM. Finally, (WIMPs) are 
considered to be decoupled while non-relativistic with a Maxwell-Boltzmann distribution function. 

Eq.JTj]) clearly shows that (WIMPs) with m ~ 100 GeV that decoupled kinetically in LTE at T - lOMeVlH 
dominate all other contributions to kf s resulting in an extremely small free streaming length A/ s ~ 6.5 x 10~ 3 pc 
unless their fractional abundance v W i„ ip < 10~ 12 . 

If (DM) is dominated by sterile neutrinos (produced either resonantly or non-resonantly) that decoupled near 
the QCD scaleJH, |33[ with m ~ keV, we find that the typical free-streaming lengths today are A/ s ~ 2 — 7kpc, 
where the larger value corresponds to the non-resonant and the lower value to the resonant production mechanisms 
respectively. The lower values are consistent with the recent numerical study of the formation of the first stars in 
filamentary structures [l^] and the values for the "cores" extracted from the data in ref.[l4| for dwarf spheroidal 
galaxies (dSphs) (r c ~ 0.5 kpc). The upper values agree with the cores extracted from the data in ref.[l5[ for spiral 
galaxies (r c ~ 10 kpc). 



II. FREE STREAMING LENGTH FOR MULTICOMPONENT DARK MATTER 



We study DM particles that decoupled in or out of equilibrium while relativistic or non-relativistic, but that are non- 
relativistic today. As argued above, Gilbert's detailed numerical study [251] confirmed by Bond and Szalay[26] shows 
that the value of k f s can be extracted from the marginal case between modes that grow under gravitational instability 
and those that are Landau damped. For this marginal case linear perturbations are stationary, and just as in the case 
of the Jeans length this marginal value can be reliably extracted in a non-expanding cosmology, including a posteriori 
the scale factor dependence of the various quantities in the Jeans wavelength which separates the gravitationally 
stable and unstable modes [Hj]. In section (|III Aj) we present arguments that determine the redshift dependence of the 
free streaming length under a suitable approximation. 

Consider several species of DM candidates that are non-relativistic today with masses m a and distribution functions 
f a where the label a = 1,2- •• refers to the different components. Each component (a) obeys the collisionless 
Boltzmann-Vlasov equation [2(| [2ll. [23|. 



df a (x,p;t) _p 

TO, 



d . _ V £ f a (x,p]t)-m a Vg$($]t)S/ p f a (x,p]t) = (2.1) 



where <I> is the total Newtonian potential which is the solution of the Poisson equation 

V 2 $(f;i)=47rG^TO Q . 9a J AjLf a ($,p-t), (2.2) 

where g a is the number of internal degrees of freedom. 

We note that whereas each individual species obeys its own collisionless Boltzmann-Vlasov equation, the Newtonian 
gravitational potential is determined by all the components as indicated by the Poisson equation (|2.2p . Therefore, 
although the different DM components do not interact directly they interact indirectly via the self-consistent gravi- 
tational potential since all of the DM components act as source of this potential which enters in Boltzmann-Vlasov 
equation of each component. 

Linearizing the Boltzmann-Vlasov equation and writing 

f a (x,p;t) = f^(p) + f^(x 1 p;t) ■ ^(x;t) = ^°(x;t) + ^(x;t) (2.3) 

where fa°\p) are the distribution functions of the species that decoupled in or out of LTE, the only assumption is 
that these are isotropic, namely only depend on p = \p\. In the non-expanding case the zeroth order equation requires 
to invoke the usual "Jeans swindle" (see the textbooks 23]) whereas in the expanding case the zeroth order equation 
is solved in terms of the inhomogeneous gravitational potential which yields the expanding background ( see [201]'). 
It is convenient to perform a spatial Fourier transform of the perturbations in a volume V 

fP{S,p;t) = ^Ff'ftp;^'^ ; $«(f;t) = -^^e*V 1 )(fc;t) (2.4) 

k k 

in terms of which the Vlasov and Poisson equations for the perturbations become 

dFi l \k,p',t) . > . n +\ d f« 0) (P) n m k\ 

\-t F a >{k,p',t) -im a k-p <pi{k;t) — =0, (2.5) 

at m a dp 
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Inserting eqn. (|2.6[) into eqn. (|2.5[) makes explicit that all the DM components are actually interacting through the 
self-consistent gravitational perturbations which act as a dynamical self-consistent mean field, as discussed above. 

Decaying or growing perturbations must be treated as an initial value problem, and following the treatment of 
Landau damping in plasmas [22| we introduce the Laplace transform of the perturbations 



(fc, p. s) = / e -*t F p (k, p; t) dt ; (k;s) = / e~ s (k ; t) dt 



The Laplace transform of the Bolztmann-Vlasov equation (|2.5|) leads to 

Z0-)ru - \ ^ ^qfc ' g #i 0) (g) . F^\k,p;t = 0) 



S + i—^- 



(2.7) 



(2.8) 



Taking the Laplace transform of (|2.6[) . multiplying (|2.8|) by (—4:irGm a g a /k 2 ), summing over a and integrating in p 
we are led to 



i47rG 



E TOa 3a 



fc 2 s) 

where the gravitational "polarization" function is given by 

,, . 4ttG V— * 9 

E(fe;s) = 1 + -p- 2^ m a5a 



rf 3 p F a w (fc,p;f = 0) 
(2tt) 3 if 



(2.9) 



d 3 P k-p dfr( P ) 

(2ir) 3 kj[ _ a dp 



(2.10) 



The collective excitations of the collisionless self-gravitating system correspond to the poles of tp^ (k ; s) in the complex 
s-plane, these are the zeroes of e(k;s). The time dependence of the perturbation of the Newtonian gravitational 
potential is obtained by the inverse Laplace transform 



ds 
2~7ri 



(2.11) 



where C stands for the Bromwich contour parallel to the imaginary axis and to the right of all the singularities of 
tp^ik] s) in the complex s-plane. 

Using that fa°\p) is a function of p = \p\ it is convenient to write 



k-p m a 



k-p 



k-p 



(2.12) 



which allows to extract the s = contribution. The resulting expression for e(k; s) becomes 



e(k; S )= £ (k;0) + iS ^^g am lJ^ J ?_ 



dfj 0) (p) 
dp 



where 



e(fc;0) = l + ^ 



J29am 3 a p 
Jo 



dti°\p) dp 
dp 2?r 2 
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In eqn. (|2.14p we have integrated by parts and introduced the free streaming momentum kj s given by 



\y2 la 



(2.13) 



(2.14) 



(2.15) 
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where 



and 



p(°>=m Q ni°> ; n g» = 9a / ^/<°>(p) (2.16) 



/J_\ f^<f(o)( v) (217) 

Kyi/a- n (°)J (27T)3 p»'» (2.17) 

the free streaming length is obtained as 



2tt 

A /s = — . (2.18) 



The expression for kf s can be written as 

<W£y /f\ Fo = 122 . 5 ^ (2 . 19) 
/s (Mpc) 2 ^ \y2/ a 

where the partial fractions are defined as 



S 



DM 



v a = 1 . (2.20) 



The collective modes correspond to the zeroes of the "gravitational polarization" function e(fc; s) eqn. (|2.13jl in the 
complex s— plane. These yield the time evolution for the gravitational perturbations 

if^(k;t) ocJV'W* (2-21) 
p 

where s p (k) are the zeroes of e(fc; s). 

It is illuminating to compare the expression (|2.15j) with that for the Jeans wavevector for a single fluid 

*3 = ^, (2.22) 

where c s is the adiabatic speed of sound of the fluid. We see that for a single collisionless component we can obtain 
the free streaming length from the Jeans length by the replacement 



\y2 



(2.23) 



which in general is different from replacing c s by the velocity dispersion y (V 2 ). This difference becomes important 
when the unperturbed distribution function favors small values of the momentum. This observation will become 
crucial when we study Bosonic particles and sterile neutrinos decoupled when relativistic but out of LTE. As it will 
be seen in detail below, in these cases the distribution function favors the region of small momentum which leads to 
dramatic consequences in the difference between 1/(V 2 ) and (1/V 2 ). 



A. Landau damping and Jeans instability 

It is clear from (I2.13[) that there is a pole in the Laplace transform of the Newtonian perturbation for the marginal 
value 

s = ; k = k fs . (2.24) 

This is akin to the marginal value k = kj in a fluid where kj is the Jeans wave vector, in a fluid for k > kj 
pressure gradients hinder gravitational collapse and the perturbations are simple acoustic oscillations, for k < kj 
pressure gradients cannot prevent the collapse and the self-gravitational fluid undergoes the Jeans instability towards 
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gravitational collapse. We can study the dynamics of collective excitations in region k w kf s searching for zeroes in 

e(k; s) for s w 0. The second term in (|2.13p can be evaluated by performing the angular integral and using that /i 0) 
only depends onp= \p\, namely 



d(cos 6) m a 
hp 



kp cos 9 
-1 — IS 



hi 



kp 



IS 



(2.25) 



For s w the branch is defined by the prescription determined by the Bromwich contour Re(s) > 0, which is 
precisely recognized as the Landau prescription for the evaluation of the integrals [22j. We find 



e{k;s) = 1 - 



k 2 

k 2 



G 

irk 3 



9a m A a 



dfi°\p) 
dp 



In 



kp 



kp 



IS 



ITT > dp 



(2.26) 



Because for small s the logarithm in (|2.26[) is linear in s, the last term in the bracket (— in) contributes to the 
leading order in k 2 



k 2 s . To lowest order in s for k w fc/ s , the condition e(k; s) = yields 



s(fc) = C [fc^ s - fc 2 ] ; C 



GLs»^/i 0) (o) 

From the time evolution of the gravitational perturbation eqn. (|2.2ip we find 



>0. 



(2.27) 



s(k) < for k > k fs 
s(k) > for k < k fs 



Landau damping 
Jeans instability . 



(2.28) 



The long wavelength limit k — > is obtained by expanding — is) 1 in the integrand in cqn. (|2.13p in powers of 
k ■ p/ras. In the resulting expression only the odd powers survive the angular integration. Keeping up to (k ■ p/ms) 3 
the long wavelength limit of e(k; s) is found to be 



a L 

and the zeroes of e(k] s) in the long-wavelength limit are found to be 

s±(k) = ± 



n 2 j-v^k 2 



where 



(2.29) 



(2.30) 



(2.31) 



and v a are the partial fractions. The Jeans frequency f2j is the same as that for single component fluids, however 
the relationship between the free streaming wavevector kf s (|2.15|) and the Jeans frequency £lj is different from that 
of the Jeans wavevector and the Jeans frequency in a single component fluid. For a single collisionless component we 
find 



i) 

V 2 ' 



kfs 



(2.32) 



whereas for a single fluid one finds [2C 



c s k. 



(2.33) 



B. An example: the Maxwell-Boltzmann distribution 



In general the momentum integral in (|2.10p cannot be found in closed form without approximating the distribution 
function. However the Maxwell-Boltzmann distribution provides an example for which a closed form expression for 
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(I2.10p can be found. This distribution function is relevant for the description of WIMPs which are heavy relics that 
decoupled in LTE while non-relativistic[l8|. In this case 



(2.34) 



where the normalization J\f is obtained from the solution of the kinetic equation for the distribution function, which 
can be found in section (5.2) in ref.[I3|. The particle density is 



P = mg 



d 3 p 
(2tt) 3 



fW{0)=Nmg 



mT d 



2tt 



(2.35) 



The momentum integrals in eqn. (|2.10j) can be carried out straightforwardly because f^°'(p) is a function of p 2 . This 
is achieved by splitting the vector p into components parallel and perpendicular to k. With d 3 p = dp^ d 2 p± the 
integrals along the parallel and perpendicular directions can be done straightforwardly. We obtain the result 



k 2 

e{k- s) = 1 - -£ + — gAfm 4 e 5 ' 
k z k a 



1 = / e _t dt 

V n Jo 



5 = 



m 
2Td 



(2.36) 



It is straightforward to analytically continue this function to Re(s) < according to the Landau prescription [22 
With /(°)(0) =N we find for the marginal case s — the solution k 2 = k 2 s , where 



k 2 s = 4nGp 



V 2 



(2.37) 



and p is given by eqn. ()2.35|) . Furthermore it is a simple exercise to confirm that near the marginal case the pole in 
the function e(k; s) is given by eqn. (|2.27ll . For the Maxwell-Boltzmann distribution it follows that 



(V 2 



3T / 1 

m ' \V 2 



m 1 



T 



(2.38) 



The long- wavelength limit can be obtained by expanding the integral in (|2.36[) for 5 1, namely 



"* dt 



1 - 



1 3 



(2.39) 



which after tedious but straightforward algebra leads to the expressions for the poles given by eqn. (|2.31|) with v = 1 
and Q 2 j = inGp where p is given by ([2~35]) and (V 2 ) given by ([2~38]) . 



III. FREE STREAMING LENGTHS FOR DECOUPLED PARTICLES 

The distribution function of decoupled particles in a homogeneous and isotropic cosmological background in absence 
of gravitational perturbations is constant along geodesies and obey the Liouville or collisionless Boltzmann equation 
in terms of an affine parameter A[H, [34], [36[ 



dX 



f[P f ;t]=0=> 



df[Pf,t] 
dt 







(3.1) 



where Pf — p c /a(t) is the physical momentum, and p c the time independent comoving momentum. Taking Pf as an 
independent variable this equation leads to the familiar form[l8|, [34], [36|] 



dt 



8P 



I 



(3.2) 



where H = a/a is the Hubble parameter and P/a = p c = constant is a characteristic of the equation. Obviously a 
solution of this equation is 



f[Pf,t] = f d [a(t)P f ]=f d \p c ] 



(3.3) 
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If a particle of mass m has been in LTE but it decoupled from the plasma with decoupling temperature Td its 
distribution function is 

f d { Pc ) = — =J (3.4) 

e T * ±1 

for Fermions (+) or Bosons (— ) respectively allowing for a chemical potential. 

Since the distribution function is dimensionless, without loss of generality we can always write for a decoupled 
particle [13] 

fdipc) = Id f |r! ( 3 - 5 ) 

where a, are a collection of dimensionless constants determined by the microphysics, for example dimensionless 
couplings or ratios between Td and particle physics scales or in equilibrium fid/Td etc. To simplify notation in what 
follows we will not include explicitly the set of dimensionless constants m/T^; aj, etc, in the argument of fd, but these 
are implicit in generic distribution functions. If the particle decouples when it is still relativistic m/Td — > 0. 



It is convenient to introduce the dimensionless ratios 34 



V= T d = TS) 5 Ut) = W) ; Xd =T d (3 - 6) 



We emphasize that the distribution functions (|3 . 5[) are general and not necessarily describing particles decoupled 
while in local thermal equilibrium. When the particle becomes non-relativistic, its contribution to the energy density 



p = m n{t) (3-7) 

where [34] 

n(t)=g^M r y 2 fd(y)dy (3.8) 

27r In 



and g is the number of internal degrees of freedom. From entropy conservation [18|, |36|, the decoupling temperature 
at redshift z is related to the temperature of the CMB today by 

T d {z) = (l + z) C^j 3 T cmb = 3 2.348 x 1(T 4 (1 + z) eV , (3.9) 

where gd is the number of effective ultrarelativistic degrees of freedom at decoupling. For a given species (a) of 
particles with g a internal degrees of freedom that decouples when the effective number of ultrarelativistic degrees of 
freedom is gd, a , the relic abundance today (z — 0) is given by[34j 



If this decoupled species contributes a fraction v a to dark matter, with £l a = v a ^-DM and taking SloMh 2 = 0.105 37 
for non-baryonic dark matter, we fmdpjij 

" a = (-^fev) ft I y2fdAy)dy - (3 - 11} 

The constraint on the fractional abundance < v a < 1 yields a bound for the mass of the particle [34j. 

We are now in position to establish a relation between the distribution function of the decoupled particles as dark 
matter candidates and the results from the Boltzmann-Vlasov equation which determine the free streaming length. 

The analysis of the Boltzmann-Vlasov equation in the non-expanding case applies to the description of perturbations 
at very low redshift in the Newtonian approximation. Such approximation is correct provided the relevant length 
scales, for example the free streaming length, are much smaller than the Hubble radius. Such is the case whenever 
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the velocity dispersion of the particles is (V/c) 2 <C 1. The momentum that enters in the Boltzmann-Vlasov equation 
is the physical momentum, which in the non-relativistic limit is related to the velocity as 



p2 

V 2 = -L 

9 



(3.12) 



The unperturbed distribution functions for decoupled particles that enter in the linearized Boltzmann-Vlasov eqn. 
(|2.5p are the solutions of the unperturbed collisionlcss Boltzmann equation, namely 



/<°>(p) =/<%,) 

where fd(p c ) are given by eqn. (|3.5[) . Restoring the speed of light c, the average 



la- 



\V 2 'a \T d:a {z) C/ 

where Td, a {z) is given by eqn. (|3.9p for the species a, and we have introduced 



Ia = 



f °°y 2 C(y)dy 



(3.13) 



(3.14) 



(3.15) 



This dimensionless ratio of integrals only depends on the ratios m a /Td, fJ-a/Td and dimcnsionless couplings from the 
microphysics at decoupling. Using eqn. (|3.9p we obtain 



(Yl\ _ L905 ( ma \ 2 J T 
\V2/ a ~ (1 + z) 2 \eVJ 9d > a *' 



(3.16) 



Inserting this result into the expression (|2.19[) and using QjjMh 2 — 0. 105 [33| for non-baryonic dark matter we find 
today, for z — 



k 2 



0.4471 2 



kpc J 



(3.17) 



For light particles that decouple while they are ultrarelativistic the distribution function fd, a {y) does not depend 
on the ratio m a /Td, however, for particles that decouple when they are non-relativistic, their distribution function is 
typically a Maxwell-Boltzmann distribution which does depend on this ratio. 

Inserting the result p. lip for the fractions p. lip and (|3.16l) into eqn. (|3.17p leads to the alternative form 



k 2 

K fs 



r5.632i 2 



kpc 



\ - 9a / m a y 
^ J VkeV/ 

a 9d,a 



fd°l(y)dy 



(3.18) 



where g a are the internal degrees of freedom of the particle of species (a) . The simple expressions (|3.17I3.18[) are some 
of the main results of this article. 



A. Redshift dependence 

The Boltzmann-Vlasov equation in a non-expanding cosmology [|2.ip is obtained by setting the scale factor a = 1 in 
the same equation in an expanding cosmology (see ref . [20l [2l| ) . Thus the linearized equation describes the properties 
of the collective excitations today. In an expanding cosmology the decay or growth of perturbations is no longer 
exponential but typically a power of the scale factor (20l l2l|. However, here we are not concerned directly with the 
manner in which linear perturbations grow or decay, but with the marginal wave- vector kf s that determines the 
crossover of behavior from growth to damping of collective excitations. If we assume that the expansion is sufficiently 
slow that it can be treated adiabatically we can obtain the redshift behavior of the free-streaming length A/ S (z) by 
replacing the densities, velocities and wavevectors with the corresponding scale factors given by eqns. (|3.8l3.9l3.16p . 
namely: 

p -> p(z) = p(0)(l + zf 




(3.19) 
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where 2 = refers to today. Defining the comoving free-streaming wavevector kt s (z) 
the corresponding scale factor 

k fs -^k fs (z)(l + z). 
Assuming the validity of this adiabatic scaling in eqn (|2.15p we obtain, 



2n/\f s (z) upon rescaling by 



(3.20) 



k) s (z) 



AttG 



V/sV ~ y (1 + z) 2 v ' ^V 2 (z) 



(3.21) 



Since the DM density and velocity dispersions for all components scale as in eqn. (|3.19p we find the following redshift 
dependence of the comoving free-streaming wavevector 



kUz) 



AttG 



E^ (0) (o) 



i 



V 2 (0) 



(3.22) 



This result is similar to the expression for the Jeans's wavevector in a Newtonian fluid[20j, |2l| upon replacing 
(l/V 2 {z)) — > c 2 {z) where c s (z) is the (adiabatic) speed of sound in the medium as a function of redshift. In fact the 
validity of this assumption is confirmed not only by the similarity with the familiar Jeans' result for Newtonian fluids 
in an expanding cosmology, but also by the exact solution obtained in ref. |29j for the case of neutrinos decoupled in 
LTE. Therefore we identify kf s given by eqn. (|3 . 1 T|) as the comoving free streaming wave- vector. The scaling behavior 
of the comoving free streaming length A/ S (z) = A/ s (0)a/ (1 + z) leads to the free-streaming mass 



Mfs(z) 



4-7T 



A/ S (z) 
l + z 



M /a (0)(l + *)* 



(3.23) 



a relation similar to the that of the Jeans' mass in the non-relativistic regime. Therefore, under the validity of the 
adiabatic assumption, the simple re-scaling of the free-streaming wave- vector and length given by eqn. (|3.22p indicates 
that we can simply obtain these quantities today (z = 0) and extrapolate to an arbitrary redshift z via eqn. (|3.22|) 
provided the redshift is still small enough that the species are non-relativistic. 

The validity of the adiabatic assumption relies on the fact that in the non-relativistic regime with (V 2 /c 2 ) < 1, the 
free-streaming length is much smaller than the Hubble radius, which is found below to be a consistent assumption, 
or alternatively kf s /H 3> 1. And as mentioned above the result for the free streaming length obtained from this 
adiabatic hypothesis is similar to the usual result for the Jeans' length[20j, |2l| and is confirmed in ref.[29( for the case 
of a neutrino thermal relic. 

A more detailed analysis of Gilbert's equation for mixtures of DM components with arbitrary (but isotropic) 
distribution functions in the adiabatic approximation will be provided elsewhere [38|. 

We now gather the above results to give the general expression for the free streaming length of an arbitrary mixture 
of non-relativistic species that decoupled in or out of LTE either ultrarelativistic or non-relativistic, in terms of the 
partial fraction (y) that each contributes to the (DM) content and the dimensionless ratios I a 



1 



(1 



0.071v 



kpc J 



E j 



| / m a y 



where 



(3.24) 



(3.25) 



in terms of the general distribution functions ()3.5() which only depend on the ratios m a /Td, Ha/Td and dimensionless 
couplings and are completely determined by the microphysics at decoupling. 

We now proceed to obtain the contributions to the free streaming wave-vector today (z — 0) from the various 
components: thermal relics that decoupled either relativistic or non-relativistic in LTE and non-thermal relics that 
decoupled while relativistic but out of LTE. 



12 

B. Thermal relics 

Let us consider ultrarelativistic Fermionic or Bosonic particles decoupled in LTE with chemical potentials and with 

a/T d , a « 1. 

• Ultrarelativistic Fermions: Neglecting m/T d in the ultrarelativistic limit, but keeping the chemical potential 
H, the distribution function is 

^ = ^Tl ; u =r d (3 - 26) 

where we have neglected m/Td <C 1. For this distribution 

/ ti°\y)dy = ln[l + e "] (3.27) 
Jo 

Combining this result with eqn. (|3 . 1 8[1 we note that larger chemical potentials lead to shorter free streaming 
scales. 

Denote If[u] the ratio I a , eqn. (|3 . 1 5|) for an ultrarelativistic Fermionic thermal relic with chemical potential /x. 
It is depicted in fig. ([1]) as a function of u = fi/T^. For fi = we find 




(3.28) 



-1.0 -0.5 0.0 0.5 1.0 



U 

FIG. 1: I F [u] vs u = y,/T d for f d {y) = l/(e (y - u) + 1). 7 P [0] = 0.3844. 

For /i = and taking v = 1 for this Fermionic species and vanishing abundance for all others, the result for kf s 
given by eqn. (|3.17[l with l a = 2 ln(2)/3£(3) agrees with that found in ref.[2!|. Therefore, a Fermionic thermal 
relic as a unique DM component yields a free-streaming length today (z = 0) 

M0) = i 1 ! (— ) kpc . (3.29) 

For example a neutrino with m ~ eV decoupling at Td ~ 1 MeV, when g d ~ 10, yields a free streaming length 
today 

A/ s (0) ~ lOMpc, (3.30) 
which is the usual estimate for free-streaming lengths for HDM candidates. 



• Ultrarelativistic Bosons: The distribution function is 
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where we have neglected m/Td <C 1 and in this case the chemical potential n < so that the distribution 
function (a probability) is manifestly positive definite. For this distribution 



roc 

/ /f (y)dy = -\n[l 
Jo 



(3.32) 



This expression clearly reveals that for m = 0, /j, d = the distribution function diverges logarithmically at 
p = 0. In an expanding cosmology with a particle horizon the smallest wavevector that can describe causal 
microscopic physics is determined by the Hubble scale. Since we are considering decoupling in LTE, it is 
consistent to assume an infrared momentum cutoff of order Hd ~ ^fg^T^/Mpi, the Hubble parameter at the 
time of decoupling. Therefore assuming such an infrared cutoff, keeping m/Td <C 1 but non-zero and with 
vanishing chemical potential we find 



ff\v)dy 



hi 



H,, 



2T d 



(3.33) 



In the denominator in eqn. (|3.15[) we can set Xd = in the ultrarelativistic limit and for vanishing chemical 
potential we find 



Ib[xcL,0] 



1 

2C(3) 



In 



2T d 



(3.34) 



Keeping the mass and the chemical potential the distribution function becomes 



for which we find 



,0), v 1 m \i 

U \V)= r~^~, — 2 ' x d = 7f- ; u d = — (3.35) 



/ f d 0) (y)dy = x d Y / e- l ^K 1 [lx d ], (3.36) 

Where K\ is a Hankel function. For m/Td, fi/Td <C 1 we can neglect the mass and chemical potential in the 
denominator in eqn. Q3.15P and we find for an ultrarelativistic boson with a non-vanishing chemical potential 
in the limit m/Td, /J,/Td <C 1 

oo 

I B [ Xd ;u d ]~^--J2 e ~ l]UdlK l[ lx -i\ ■ ( 3 - 37 ) 
AW l=1 

This function features an infrared divergence in the limit m/Td, fi/Td — * from the numerator in (|3 . 1 5|) given 
by (|3.36j) . which must be regulated by assuming an infrared cutoff of the order of the Hubble scale at decoupling 
as in eqn. (|3.33p . A more thorough analysis of the infrared behavior of the decoupled distribution function 
in an expanding cosmology is necessary, in particular an assessment of the (causal) thermalization process for 
super horizon wavelengths. Such program, interesting all by itself, is clearly beyond the realm of this study. 

Bose-Einstein Condensation: If chemical freeze out occurs before kinetic decoupling, it is possible that bosonic 
particles undergo Bose-Einstein condensation [34l| . Under this circumstance the homogeneous condensate (p — 0) 
contributes to the dynamics of the scale factor and inhomogeneous perturbations only modify the distribution 
function of the particles outside of the condensate. The linearized Boltzmann-Vlasov equation therefore applies 
to the non-condensate part, since the dynamics of the condensate cannot be described by the incoherent Boltz- 
mann equation but by the equation of motion of the coherent and homogeneous condensate, which in turn is 
coupled to the Friedmann equations for the scale factor. 

In the case of Bose-Einstein condensation the chemical potential attains its maximum value fj, — m in which case 
the contribution from the particles outside the condensate is lB[xd,%d] which has a strong infrared divergence 
for any value of the mass. This is a consequence of the fact that for ultrarelativistic Bose- Condensed particles 
the distribution function diverges at p = for any value of the mass. 

However, as discussed above this infrared divergence associated with a Bose-Einstein condensate must be care- 
fully assessed in an expanding cosmology because infrared modes with wavelengths larger than the size of the 
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horizon will not be in LTE since no causal processes can establish thermalization for superhorizon wavelengths. 
Hence we conjecture that the integrals of the Bose-Einstein distribution function must be cutoff in the infrared 
at a momentum of the order of the Hubble scale at decoupling, H d ~ JlfdT^/Mpi. Such a cutoff leads to the 
estimate 



Ibec 



mMpi 10 
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g2Tj Jg2 VkeW V T d 



m \ /GeV\2 

Tev) \~fr) 



(3.38) 



If this were the only DM component, the resulting comoving free-streaming length at z — (today) is given by 



A/.(0) 



0.014 pc /keV\f / T d \ 

„& \~) \GeVJ 
9 d 



(3.39) 



It is clear from the discussion above that the microphysics of decoupling of light bosonic particles, with or 
without a Bose-Einstein condensate requires a thorough assessment of the infrared behavior of the distribution 
function. The primordial velocity dispersion is very sensitive to this cutoff whose origin lies in the causal aspects 
of decoupling. This important physical aspect must be studied in deeper detail, a task that is certainly beyond 
the realm of this article, but we can nevertheless conclude that a light bosonic particle decoupled in LTE can 
effectively act as CDM as a consequence of the infrared sensitivity of the moment (1/V 2 ) for Bosonic thermal 
relics either condensed or not. 

• Non-relativistic particles: The distribution function after freeze out in LTE is the Maxwell-Boltzmann 
distribution 



mT d 
2vr 



rid 



mTd 



2- 



where n d is the number of particles per comoving volume at freeze-out [l8j 

2ir 2 



n d 



45 



(3.40) 



(3.41) 



and Yqo is obtained from the solution of the kinetic equation and is a function of the annihilation cross section 
(see section 5.2 in ref.[18|). We find 



fd°\y)dy 



45 x d 



(3.42) 



and for the integral I a eqn. (|3.15[) denoted by I nr. for the non-relativistic (Maxwell-Boltzmann) distribution 
we obtain 



J-nr — — 



Td 
m 



(3.43) 



In the case of (WIMPs) with[35j, m ~ 100 GeV ; T d - lOMeV as candidates for cold dark matter, x d ~ 10~ 4 . 
If this is the only DM candidate with v = 1 with vanishing abundance of the other WDM or HDM candidates, 
the comoving free streaming length at z = is given by 



A/s(0) 



0.014 pc / 100 GeV \ \ /lOMeV 



91 



^iuuuev^3 ^ 



■X 



(3.44) 



from which it follows that for WIMPs with m 
g d - 10 [l| 



T d ) ' 

100 GeV that decouple kinetically at T d - 10MeV>5] when 



A /s (0) - 6.5 x 10~ 3 pc 



(3.45) 
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C. Decoupling out of LTE 



In ref.[34| the following distribution function for particles that decouple out of LTE and that effectively models 
several cases of cosmological relevance was introduced, 



fd(y) = /o/eg(^J %o - y) 



(3.46) 



where f e q\Sjr) is the equilibrium distribution function for a relativistic particle at an effective temperature rfTd- This 
form is motivated by detailed studies of production [3i| and thermalization process that proceeds by energy transfer 
from long to short wavelengths via a cascade with a front that moves towards the ultraviolet [40] . If the interaction 
rate for mode mixing becomes smaller than the expansion rate the advance of this front is interrupted at a fixed 
value of the momentum, identified here to be p° c = VnTd where Td is the temperature of the environmental degrees of 
freedom that are in LTE at the time of decoupling [34]. The amplitude fo and effective temperature rfTd < Td reflect 
an incomplete thermalization behind the front of the cascade and determine the average number of particles in its 
waA;e[40J. The non-equilibrium distribution function (|3.46p yields a fairly accurate description of these processes and 
the decoupling out of LTE. 

Remarkably, this non-LTE distribution function also describes [HJ sterile neutrinos produced non-resonantly via the 
Dodelson-WidrowJUJ (DW) mechanism or resonantly via a lepton-driven MSW resonance [33j. 
For the general form (|3.46[) of the distribution function we find 



r/ z 



Vc 



riT d 



H[s] 



Jo fe q {y)dy 

Io S y 2 feq(y)dy 



(3.47) 



The function H(s) for f eq (y) = l/(e y + 1) is a monotonically decreasing function of s with limiting behavior H(s) 
3/s 2 for s -> and H(s) -> 21n(2)/3C(3) for s -> oo, it is displayed in fig. © in the interval 0.25 < s < 2. 



FIG. 2: H(s) vs s for f eq (y) = l/(e y + 1) 



Although a more detailed assessment of the production and pathway towards thermalization of sterile neutrinos is 
still being debated [ZJ H3] we will adopt the semi-phenomenological results of refs.JH, [33[ as guiding models for the 
distribution functions of sterile neutrinos decoupled out of LTE. 

• Sterile neutrinos produced via the (DW) mechanism 32]: for this case rj = l;s — > oo; fo ~ 

0.043[keV/m] [32l 1341] . The integral I a does not depend on the amplitude fo (only the fractional abundance 
v is affected by fo) and is therefore given by the value for an ultrarelativistic fermion with vanishing chemical 
potential decoupled in LTE[32], namely low = 21n(2)/3C(3) ~ 0.3844. // these sterile neutrinos are the only 
DM component, namely v = 1, the free streaming length today is given by the same form as for a Fermionic 
thermal relic eqn. (|3.29p but with vanishing chemical potential, 



A/ s (0) 



22.7 /keV 



fi;l 



v m 



kpc . 



(3.48) 
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In the Dodelson-Widrow scenario[32j the sterile neutrino production rate peaks at T ~ 130 MeV which is near 
the QCD scale. Taking the decoupling temperature in this range results in gd ~ 30[l8j], leading to 

/keV\ 

A /S (0)~7kpc (—J . (3.49) 

Therefore m ~ keV sterile neutrinos produced via the (DW) mechanism yield free streaming lengths today of 
the order of A/ s ~ 7kpc. As discussed in ref.[13] there is a potential ambiguity in gd because near the QCD 
phase transition there is an abrupt change in the effective number of relativistic degrees of freedom. However, 
because the cube- root of g& enters in A /„ the free streaming length is not very sensitive to this ambiguity. 

• Lepton-driven resonantly produced sterile neutrinos [33j: for this case r\ = l,s ~ 0.7; fo = lJU, HH, and 
we find I a — H(0.7) = 6.814. If this species is the only DM component we find the free-streaming scale today 

A /s (0) = — ( )kpc. (3.50) 

The production rate in the resonant case also seems to peak near the QCD scale[33] at which gd ~ 30. Taking 
this value for the decoupling temperature we obtain the estimate 

A /s (0) = 1.73 kpc ( — -J. (3.51) 

Hence m ~ keV sterile neutrinos produced via the lepton-driven resonant mechanism [33l| mechanism yield free 
streaming lengths today of the order of A/ s ~ 2 kpc, which is consistent with the values ~ 3 kpc used in[lj| to 
study the first stars forming in filamentary structures. 

The enhancement of I a and consequently the shortening of the free streaming scale in the case of sterile neutrinos 
produced resonantly out of LTE follows from the fact that the distribution function resulting from the resonant 
production mechanism favors low momenta. Therefore for the same value of the contribution of sterile neutrinos 
produced resonantly or non-resonantly out of LTE, the free streaming length in the case of resonance production 
is ~ 4 times smaller than either a thermal Fermionic relic with vanishing chemical potential or a sterile neutrino 
produced non-resonantly via the (DW) mechanism. Just as in the (DW) scenario, there is an ambiguity in the 
precise determination of gd near the QCD scale, but again the free streaming length is not very sensitive to this 
ambiguity because of the power 1/3. 



Although the details of the mechanism of production and decoupling of sterile neutrinos are still under scrutiny 41| , 
we take the above results as a guideline. In particular the lesson from the resonant production case highlights the 
potentially dramatic reduction of the free streaming length when the non-equilibrium distribution function favors 
smaller values of the momentum. This observation clearly calls for a deeper assessment of the production and 
thermalization of sterile neutrinos to extract reliably their non-equilibrium distribution function. 

Combining all the results above into eqns. (|3.24I3.25[) . the free streaming length for a mixture of Fermionic and 
Bosonic thermal relics, including a possible BEC, WIMPs and sterile neutrinos produced non-resonantly or resonantly 
assumed to be described by the respective distribution functions quoted in refs.[13, [H[ is given by 



1 1 



y fs (z) (1 + z) 



0.071i 2 

species 



kpc 



\ 1 VF 9lp V keV ) IF[U] + Vs 9 *> s V keV ) 6-814 + VB 9 ^ B \ keV ) Ud] + 

10 12 u wimp gl mrnp ( 10Q ^y) (tTT^f) I ' ( 3 ' 52 ) 



where the index F refers to thermal Fermions and sterile neutrinos produced non-resonantly via the (DW) 
mechanism 32] for which If[0] = 21n(2)/3£(3) = 0.3844, and the label s refers solely to sterile neutrinos produced 
resonantly via the mechanism in ref.[33|]. In each case g^ a is the effective number of ultrarelativistic degrees of freedom 
when the corresponding species decouples. For sterile neutrinos produced by cither mechanism gd ~ 30 corresponding 
to a decoupling temperature near the QCD scale. 
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IV. CONCLUSIONS AND FURTHER QUESTIONS 



In this article we have implemented methods from the theory of multicomponent plasmas to study free streaming 
of a mixture of non-relativistic DM candidates that include Fermionic and Bosonic particles that decouple in LTE 
while relativistic, including the possibility of a Bose-Einstein condensate, heavy thermal relics that decoupled in 
when non-relativistic (WIMPs), and sterile neutrinos that decouple out of LTE when they are relativistic. We solve 
exactly the Boltzmann-Vlasov equation for the gravitational perturbations in a non-expanding cosmology and obtain 
the "gravitational polarization function" whose zeroes determine the dispersion relations of the collective excitations 
of the self-gravitating collisionless gas of particles. The free-streaming wave vector is obtained from the marginal 
solution that separates Landau damped short wavelength perturbations from unstable collective modes. We obtain 
the free-streaming length for arbitrary (but isotropic) distributions of the particles that decoupled in or out of LTE 
solely in terms of the fractional abundance of the different species and integrals of their distribution functions which 
depend on the microphysics at decoupling. 

Because all of the components are non-relativistic and assuming that the expansion is slow we provide an adiabaticity 
argument that allows us to implement a simple rescaling of densities, velocities dispersion and wavelengths to extract 
the redshift dependence of the free-streaming length. The validity of this adiabatic approximation is confirmed by 
the similarity of the result to the Jeans' length for Newtonian perturbations in an ex pan ding cosmology and by the 
explicit result for the free-streaming wavevector for thermal neutrinos obtained in ref. [291 ] . 

The main result for the free streaming length as a function of redshift for an arbitrary mixture of DM components 

is 



1 



1 



(1 + *) 



0.071V 



kpc J 



y * 9d.a 



( m a y 
VkeV/ 



la, 



where 



f °° v 2 &v)dy 



(4.1) 



(4.2) 



is a ratio of integrals of the distribution functions that only depends on the m a /Td, Ha/Td and dimensionless couplings 
and is completely determined by the microphysics at decoupling. Evaluating these integrals for thermal Fermionic and 
Bosonic relics (with or without condensation), WIMPs and sterile neutrinos decoupled out of LTE either resonantly 
or non-resonantly with the distribution functions obtained in refs.(32l. |33| respectively, we find the general result 



^fs(z) (1 + 2) 



0.0711 2 



kpc J 



\ I VF 9 ^ F VkeV ) IF[U] + Vs ^ s VkeV ) 6 ' 8W + VB g ' lB VkeV ) Ib ^ Ud] 

species ^ 

9d,wim P VioOGcV/ V 



10 12 v u 



.lOMeW 



where v a is the partial fraction of each component with ^ Q v a = 1, the sum over Fermionic species (F) includes thermal 
relics and sterile neutrinos produced non-resonantly via the (DW) mechanism [32| . for which the chemical potential 
vanishes and Lf\0] = 0.3844, the label s is for sterile neutrinos produced via a lepton-driven (MSW) resonance as 
described in ref. [33], and the label B stands for condensed or non-condensed Bosonic thermal relics. 
This expression features several important consequences relevant for large scale structure formation: 

• A non-negligible v wimp ^> 10~ 12 , for a CDM candidate with m ~ 100 GeV which decoupled at T<j ~ 10 MeV(35j 
overwhelms all other components (but for a possible BEC) and leads to small free streaming lengths consistent 
with CDM regardless of the presence of WDM or HDM components. For v W i mp 3> 10~ 12 the free-streaming 



length of mixed DM is completely dominated by WIMPs and is given today by 

10MeV\ 3 



A/ s (0) 



0.014 pc /lOOGcV 

\ TO 



9d 



T d > 



(4.4) 



from which it follows that for WIMPs with m ~ 100 GeV that decouple kinetically at Td ~ 10 MeV[35( when 
9d - lOpl 

A /s (0) - 6.5 x 10- 3 pc. (4.5) 

This cut-off scale might well be related to the smallest non-linear structures found in[l(| unless there is some 
substantial violent relaxation and merging. 
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• For vanishing chemical potential Ud = 0, non-BEC Bosonic ultrarelativistic relics feature an infrared enhance- 
ment in Ib [xd, 0] which must be regulated by assuming that the integrals of the distribution function are cutoff of 
order of the Hubble scale at decoupling Hd ~ y/gdT^/Mpi. If this is the only DM component, the free-streaming 
length today is 



A/ s (0) 
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kpc ; I B [x d ,0] 
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2T d 



H,, 



VW 



(4.6) 



For BEC Bosons Xd = Ud the function lE,[xd,x d ] is divergent as a consequence of the infrared divergence in 
the numerator of I a (eq. 13. 1 5j) . However, we highlighted that the physics of BEC formation in an expanding 
cosmology must be assessed in greater detail in order to understand the behavior of superhorizon modes. This 
observation also holds for the non-condensed Bose gas decoupling when relativistic with fi/Td <C 1 because the 
distribution function in this case is also infrared sensitive. In both these cases an estimate of the free streaming 
length may be obtained by introducing an infrared cutoff in the integral in the numerator of I a of the order of the 
Hubble scale, since superhorizon modes cannot establish thermal equilibrium via causal processes as discussed 
above. If a BEC is the only DM component, its free-streaming length today is approximately given by 
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Therefore even when these relics decoupled when they were ultrarelativistic, they could effectively act as CDM 
components. This is a consequence of the fact that the distribution functions favor small momenta. 

If sterile neutrinos that decouple out of LTE near the QCD scale produced either non-resonantly via the (DW) 32] 
mechanism or via a lepton-driven MSW resonance[33j near the QCD scale are the only DM components we find 
the following free-streaming lengths today 



A/ s (0) ~ 7kpc ( 1 hod — resonant 
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(4.9) 



The smaller values of the free streaming length are compatible with those in ref . [19( where it is found that first 
stars form in filamentary structures with length scales of the order of the free streaming scale and within a 
factor ~ 3 — 4 seem to be also consistent with the "cores" resulting from the fit of the density profile for dwarf 

0.5 kpc. The larger values are consistent with those found in ref.[l5[ for the 
-5- 10 kpc. 



spheroidal galaxies in ref.[14( r. 
cored profiles of spiral galaxies 



We believe that these results lead to a significant advance in the understanding of collisionless (DM) because 
trying to obtain a reliable estimate of the free-streaming lengths via the numerical integration of Gilbert's equations 
for a combination of arbitrary distribution functions corresponding to particles that decoupled in or out of LTE is 
undoubtedly a daunting task. 

The sensitivity of the free-streaming scale to the details of the distribution function at low momentum and the 
importance of a reliable determination of the free-streaming length as a measure of the cutoff of the power spectrum of 
linearized cosmological perturbations require a fundamentally sound understanding of the microphysics of production 
and decoupling of sterile neutrinos, a program currently underway [HJ. 
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